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ABSTRACT 

Trie commonly-used " 4 / 2 -Earth Model ^ in determining 
radar altitudeSj has recently teen shown to over-estimate 
curvature at elevations greater than several kilometers. 

It is further limited in that provisions are not made to 
account for commonly-occurring super-refractive atmos- 
pheric layers. Based on these shortcomings, a new model 
is constructed which incorporates an exponential refractiv- 
ity model proposed by recent investigators, together with 
superimposed low-level trapping layers of various inten- 
sities. From tils model, tables of heights are computed 
for altitude-determining radars. These heights should be 
considerably more representative than those based on the 
"4/3-Earth Model!' 

The author wishes to express his appreciation for the 
assistance of Professor F. L. Martin and Mr. E. 'Ward 
of the U . S. Naval Postgraduate School in this investi v a- 
tion. 
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SYMBOLS AND ABBREVIATIONS 



fie 

A 

/Sc*J 

h 

Ah 

n 

N 

z 

S 



Radius of the earth 
Effective radius of the earth 

Angle of ray emission, relative to the horizontal 

Ray-entrance angle, relative to the horizontal 

Ray- emerg ence angle, relative to the horizontal 

Height of trapping layer 

Thickness of trapping layer 

Refractive index of atmosphere 

Refractivity of atmosphere 

Height, measured from surface of the earth 
Rang e 

Relative curvature of ray 
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1. Introduction 



Atmospheric refraction of radar rays can seriously 
affect the accuracy of surface-hased hei T ht-findinp radars. 
The purpose of the paper is to consider some models of 
ano^olous atmospheric refraction, and determine corrections 
necessary to observed radar altitudes in order to provide 
an accurate target altitude. 



T 




r 

Figure 1. Typical radar ray path. 

Consider a transmitter, T, in Figure 1, and its target, S. 
Because of atmospheric refraction, the radar ray path will 
be, in general, a curved arc, TS. In order to determine 
accurately the altitude of the target, S, we must be able 
to predict correctly the curvature of the arc TS along all 
points of the arc. 

In 1933, Schelling, Burrows, and Ferrell C4] , set 
forth the effective earth's radius model to account for 
the effects of atmospheric refraction, better known as the 
"4/3-Earth iuodel". This is to say that if we were to depict 
the earth as a sphere of some radius suitably larger than 
that of the real earth, then normally downward tending 
radar waves could be drawn as straight lines. The accepted 
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standard effective radius has become 4/3 the radius of the 
earth. 

Recent investigations by Bauer, itason, and Vvilson^lj, 
of Lincoln Laboratories^ and Bean and Thayer [. 2D of the Rational 
Bureau of Standards, establishes that the lapse of refrac- 
tive index, dn/dz, wl ere n is the refractive index and z 
the height above the earth, is not constant, as it was 
assumed to ' e in the 4/3-earth model. These investira- 
tors suggest that the revised standard n-profile is mors 
suitably described in the mean by an exponential curve, 
rather than a linear one, such that: 

-n = 1 + C,e- C * 2 (l) 

where Cp and Cp are empirically determined constants. 

At this point, let us consider, from a Quantitative 



viewpoint. 


the 


effect 


of the ind 


ex of refract 


ion 


or a 




radar 


ray . 
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generally 


of the order 
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1.0CC3, 
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with the average or "standard" curvature, r :aultin from 
the 4./2-earth model. For-^ = 47.3 units/lOCO feet 
Equation (2) shows that the radar-ray curvature is ti c same 
as t v at of the earth. This condition is known as a trap- 
ping, or ductinp condition. Stnce vertical N-rradisnts 
well in excess of 48 units/lCQO feet exist in the atmos- 
phere, the radar-ray oath can be si riificantly alt 'red 

i* 

d urine its passage ttrourh the atmosphere. To illustrate 
the error this tendin * introduces in altitude determining 
radars, consider a transmitter T situated below a super- 
refractive layer. in Pi ure 2. 
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Fir ure 2. Example of trapping . 

The solid line TS, is the actual ray path between the 
target and the transmitter. The dashed line, TS 1 , is the 
ray path alon " which one assumes the radar energy ! as travel- 
ed, if the 4/2-earth model applies, and from which the set 
computes its displayed hei: ht. Therefore, the radar set 
will display a height IP, which, as can be seen, constitutes 
a sizeable departure from the true heipht. If. 
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L. C nutation model. 

..s explained Ir iecticr 1, the 4/3- r 
1* :r rontly used to most cp-ratlonal radons, ‘ ■ 

►vlf. respect to doth lean an., transient enrol . t .. 

•. e, ore. One cannot expect to ? j a tie tc co 
el ' ieh will account lor all th e transiet . e 
aw.ospr ere. .» model can tc constructed, ■ v . r, 
will account for certain si nil leant, rut]... r st.- ' 1. , 
p’ eric cone it lens, such as trapping. .it; this 
ind, the model depicted in ure 1 has he .n ce ■ r c • 




ri< ■ -c 3, ...cael ntn.osprere wit]', trapping Ir.r os 



An exponential lapse of n with hei. ht r s he .n sssuneu 
loth below and above the tranpin ' layer, follow in Location 
(1). A sin: le trapoin~ layer has been assumed, oh thickness 
Ah, and exhibitinp a linear lapse of n. As can be seen 
from the eometry of Fipure 3 and Eruation (1), the value 
of n at any point belo w the base, 2=k , of a trapping layer 
is expressed: 



on. = i c, e 



C,2 



o ^ 



£ < k 



(4) 



similarly, the value of n at any point within the trapping 
layer of thickness Ah can be expressed as: 

«.= (l + C,e" Czh ) + ZTT (2- h ) k4 2 «k + ak ( 5 ) 
Finally, the value of n above the trapping layer can be 
expressed: 



n = 1+ Ce' Cz * - C e' 






e 



■QaK_ | 



An. 

L\y 



+ -T- aH 



1 > K 4 A k ( 6 ) 

Note that the model ;ust constructed ’"as t^e same slope 
dn/dz above the trapmin- layer as it would have had in 
the absence of a trappin~ laj/er. V/h ether this is f lfilled 
in practice is not known; however, this feature of tie 
model is very convenient from a computational viewpoint. 

This feature may be stated anotrei- way: for any Lven 
trapping layer of specified thickness, rad lent and hei ht. 



Eouation (6) can be seen to have the form 

w = 1 v C,e" ClZ + \i(t) (7) 

where K(t) is a function of the properties of tie trap- 
inf layer, and is a constant for any set of computations 



5 



deal in 'with any sin 1? *r ppin layer. 

Levin' described t- model, the next c U, *- 
selection of appropriate values for C]_ and Co. T 1 ic.v sti 
rations of D ean and Thayer Q 2\ have established v lue _'or 
tie constants based cn 39 station-ye rs of radiosonde 
otss . vationSj reoresentin. an all— season averaf c at 13 
climatically diverse locations. Cn the other >-rj, "nuer, 
Lsson and si Ison [l] have su -ested values fc- the c'r starts 
b' sod on a selected set of soundin'S representirr a "cool, 
dry" spring atmosphere. «»hile the work of dem nr i T> cyer 
represents a r. ore precise representation of an. av^r/e "all 
cases" sounding , the att osphere described by Bauer, f-asen, 
and hi Ison better approximates an Ideal atmosphere into 
which anor.olous conditions can l e interjected. Therefore, 
the values proposed by Bauer, Lason, and 17 i Ison Lave been 



psed to refine the exponential atmosphere, and donation (1) 
now becomes: 

, -.OS'?!?. 

r\ - l"h . 0003 Z 0 ^ 

where z is .;.e: sired ir. kilo foot. 
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Ray ..some try. 







a. 



Fiyure 4. Ray eometrv. 

Consider the ray in Firure 4 enterin' •' spherically 
homogeneous layer of the atrosphcre with entrance anfle^ , 
and era e rp i nr from the top of the layer with on Is of ener- 
f ence^ +J . Annies^, and are always considered to he 
measured from the local horizontal. Snell's law tells us 

[ 3 ] : 

!/ -i (&.+ COS ^ — ^L> I h ~ L'l) COS — Conxtsnt (9) 

where n^_ and n- L+1 are the indices of refraction at tot torn 
and top of this layer respectively; and are the 
heights of the points a hove the surface of the earth; 
and a is the true radius of the earth^. lutpressin (0) 

2 For tie purposes of this paper, the earth is assumed a 
perfect sphere of radius 20SC6.9C8 kilofeet. 
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in t.rms of the logarithms of the values, ana aiffer- 



ertistin' with respect to height: 

"n if ounz Cos P > 4:a ~ ^ 

dearrang inr terras, and noticing that = ft -ft — -ft ^-3 

* ' cU r .U cU dp fiJsn p 

tie last e oust ion may be written: 



do. - 



dLjg 



( 10 ) 



CoS fi [n3ti + ^h] 

Heplacing the differentials with a finite-difference form, 
( 10 ) can be written: 



AS L 



ftn “ ft 



c 



' r*‘ |_Kti +n l X5i + i ~*i) M a + CIl 



2 - 1 



( 11 ) 



bsing Equation (S) to solve forft m -ftft one obtains: 



Cos ’ 1 



A-i -A 

Substituting (12) into (11) ^ives: 



h-L f o, -t- Cos Bj, ' 

•Vi (a.-+Z L+a ) ) 



'ft 



(12) 



AS- t = 



ft .ft (0-vlO Cos 8 i A*] D 

L Gs \ Hc , i ( ol ->-Z iv1 V JJ ft 



Cos ft 



2(n u i-aQ ^ 

( n in + ^ 0 ( 2 ; Q-* (g,-+Z,,i ) 



4- 



(13) 



^ration (13) can be numerically solved by the following 
steps : 

1 . select small increm ents A>ft — ft +1 — ft , iven an 
initial eng le ft 

2. compute ASi, and yft +1 over small thicknesses using 
Equations (4), (5) , ( 6 ) , (12) , and (13) 



3. compute range S and hei ht 2 where: 

v _tL 



s = £ ^ 



AZ t 



L --1 



(14) 



8 



^or computations used in this pap., r, A2 L r-o ' n :’x.5ir ;s 
dCC feet. 

Let us look for- a moment at bcuation (lb). if v\ ^ 
is sufficient rreater than , and the ray h -s c low 

anrle of elevation so that Cos &; 1 , t v 'e t-m 

will become reater than unity. under these conditions tie 
cosine of an anple cannot exist, and Snell's law fails. 

Such a condition represents complete reflection of t'^ -• ray, 
or trapping. Under these conditions, the Lei ht of any 
tar. et detected by this ray will be less than th • -i ht 
cf the top of the trapping layer. 
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4. Computation procedure 

Tire model as described in Section 2 was computet V 
the top of the trapping layer at every 1000 feet up to 
5000 feet, and ter thicknesses of 200, 400 and eCG feet, 
with h-rrsdients of 50, 100, and 20C units/1000 fe't. 

These values have teen selected so t v >at all co monly- 
occurrinp atmospheric trapping layers will lie within the 
ranpe of parameter values selected, providin'" a basis for 
tabular interpolation. Selection of the number of parameters 
used has been a compromise between complete tabular list in. s 
and the size of Appendix 2. The author has tried t u select 
the cases so that both reouirenents are adeouately -.at, 
keeping Appendix 2 to such a size that Its value will not 
be lost it: its bulk. 

The procedure outlined in Section 5 was pro* rammed 
for the 000-1504 hifh speed digital computer. days were 
computed for increments in the initial anfle of one- 
ei*hth decree up to one-half decree; then in hs.lf-derree 
increments from 6 to 20 decrees, then 5 decree increments 
for arifles froir 25 to 50 decrees. ITei hts were computed 
in 200 feet steps, and for each height the ranpe computed, 

T v e computed ran' e and the known initial anples were then 
used to solve the 4/5 earth height usinp the Donation Cl - }: 

^ = SS„(3 0+2 4 e (15 > 

where CL e is the effective earth's radius, in this case (4-/3)a. 
Louation (15) is notec to be an approximation formula, but 



10 






ij 









it constitutes the coi. only employed ac^.i elation ior.. fls. 
in the design of altitude-detem Lninp radars. wince t .is 
paper concerns corrections to raoar observations, rather 
than corrections to the 4/3-eartb model, the a - p oximat ion 
is appropriate. 

It will to notes that ter each simultaneous computation 
of a complete ray through the model atmosphere errs the 4/3= 
earth atmosphere, both rays have a common initial an* le and 
rartpe. The common rar e value for the two rays was induced 
to intervals of t :n miles, and the correspondin' I ts 

for each model were interpolated to proper values for the 
1C -mile intervals for identical values of . Aft t r all 
rays 1 ave teen computed, two sets of tables will exist for 
each ran e, one set containing a family of hei hts for the 
4/3-earth atmosphere, and tie second table containing corre- 
sponding values of the model beipht. The table containing 
the 4/3- earth model heipht was then reduce d to 1C, ^ CO foot 
intervals, and tr.e corresponding true hoi hts were linearly 
interpolated to their appropriate values. The p -neral 
technique can be described in general terms ss two-way 
tabular interpolation, employ inp the search and le ic functions 
of the CDC-1G04 compute r A The results were then printed in 
the form displayed in Tables 1-45. 

For esses wnere ray trappirp occurred, the computer 

^’The author has ma 5 ntained a complete set of th _'A. c?rus 
used in the main pro' ram, available t.c those uv • st s in 
tli* Precise p .* 0 : rarer. inp technicues. 
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generated ranges to the pc int of trapping ano V see ■ ire 
placed in a ranj s table opposite the 4/3- earth -i ht. 

Kei r 'ts were obtained with their coi’responding 4/3- earth 
model heights. The actual height of any trapped ray was 
arbitrarily defined as the mid-point height between the 
surface and the top of the trapping layer, ho further- 
computations were made for any ray in the model atmosphere 
after it reached the level of trapping as .just defined. 

The ray corresponding to the next increment of /? 0 was then 
computed as outlined in Sections 3 and 4. 
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5. Conclusions. 



Test runs of the computation scheme for a strict 
exponential model wore in close agreement with the publish- 
ed results of Bauer, ha sen, and .vilson, beinp' within i 5CC 
feet cf these results. The tables in Appendix I should 
therefore provide adequate accuracy for most radar uses. 

At- this point, the author would like to caution the 
reader that, while the digital computer provided results to 
eipht decimal places, at extreme ranges and lew altitudes, 
the tabl es probably do not provide accuracy to arm/ g reater 
decree than f 500 feet. In order to provide preater accuracy, 
smaller height increments and more sophisticated integration 
procedures would have to be used. 

The tabulated values in Appendix I should provide a 
valuable addition to the heigbt-f indinp" radars. The ex- 
istence of trapping layers can be ascertained usinp pro- 
cedures outlined in various AiiCV.'A publications L5]. 

This problem constituted the first computational pro- 
pram of any size completed on any of the CDC-1604 computers. 

As a result, the task was undertaken prior to the existence 
of an adequate subroutine library for the computer. Since 
the author not only had to program the primary routine, as 
well as all subroutines, an approach was selected to keep 
these subroutines to a minimum. For example, a natural loga- 
rithm subroutine and other types of interpolation subroutines 
could have been useful, but limited time precluded the develop- 
ment cf such routines durinp this particular study. 
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API JJ.A 1 



■Use of Tables 

Lpon selecting the proper table for use, based on pre- 
vious knowledge of atmospheric conditions, enter the left 
hand column with the observed ranpe, and enter across the 
top with the observed radar height, in thousands of feeto 
The true altitude of the tar-pet may then be read from the 
table . 

EXAMPLE: Assume we are usinp Table 1, with an ob- 

served heirht of 20,000 feet (20 kilofeet) and an 
observed ranre of 150 miles. The table fives : 
+.20612003 -1-02. bote that the value is listed as a 
fractional coefficient with any accompany inn power 
of ten. Movinf the decimal point two places to the 
rijit as indicated by the +02, the value becomes 
20.612003 kilofeet, or 20,612 feet. Considerinf the 
accuracy of the computations , this would e recorded 
as 20,500 feet. 

In the event trapping occurred, a notation^ has been 
included at the bottom of the table indicating the first 
anfle, of t hose angles tested , which passed threufh the 
layer without be inf trapped. Elevation angles less than 
this anfle are to be treated with caution, as the heifht 
values in the table will indicate. 



^Except if trappinf does not occur for any initial ray 
an; le ; then no notation appears under the table. 

15 
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APPENDIX 2 



TABULAR it AH CE -HE 1 3 II T VALUES 
Tables 1-45 
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HEIGHT DETERMINATION FOR MODEL ATMOSPHERE WITH TRAPPING LAYERS 
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